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We report a measurement of the exclusive e + e — » /j(*) ± D*T cross section as a function of center- 
of-mass energy near the D*-*' D* T threshold with initial state radiation. A partial reconstruction 
technique is used to increase the efficiency and to suppress background. The analysis is based on a 
data sample collected with the Belle detector at the T(45) resonance and nearby continuum with 
an integrated luminosity of 547.8 fb _1 at the KEKB asymmetric-energy e + e~ collider. 

PACS numbers: 13.66.Bc,13.87.Fh,14.40.Gx 



with charm meson pairs are of special interest because 
they provide information on the spectrum of J PC = 
1 charmonhim states above the open-charm thresh- 
old, which is poorly understood [![. The observation 
of the charmonium-like F(4260) state, seen only in 
J / ipTTir(KK) final states @, S i, H, has stimulated re- 
newed interest in this field. Curiously, the y(4260) peak 
position is close to a local minimum of the total hadronic 
cross section [gj]. The large branching fraction to J/ipTTir 
inferred from the total hadronic cross section is unex- 
pected for a conventional charmonium state with such a 
large mass and width. In a recent measurement 0] the 
e + e~ — ► DD cross section is described by known charmo- 
nium states without a significant Y"(4260) contribution. 
A study of the production cross-section of the charmed 
meson pairs in this energy range could help clarify this 
intriguing situation. 

In this paper we report a measurement of exclu- 
sive cross sections for e + e~ — > D* + D*~ and e + e~ — > 
D + D*~ [§] using initial state radiation (ISR). The data 
used for this analysis correspond to an integrated lumi- 
nosity of 547.8 ftT 1 collected with the Belle detector 
at the T(45) resonance and nearby continuum at the 
KEKB asymmetric-energy e + e~ collider [lpj ] . 



of a hard photon. The selection of e+e~ — > flW+ D*~y lsI 
signal events using full reconstruction of both the D^*> + 
and D*~ mesons, plus the 7i sr photon, suffers from low 
efficiency due to the low reconstruction efficiencies, 
small branching fractions and the low geometrical accep- 
tance for the 7isr, which tends to be emitted along the 
beam line. Here, we use a method that achieves higher 
efficiency by requiring full reconstruction of only one of 
the £>(*)+ mesons, the 7i sr , and the slow 7T^T ow from the 
other D*~ . In this case the spectrum of masses recoiling 
against the £)W + 7i sr system: 



M rcc (£ W+ 7isr ) = yJ(E C M-E DW+ ^J< 



-p 



(1) 



hadronic cross section at y/s 
3 + e~ center-of-mass (CM) en- 



peaks at the D*~ mass. Here Ej)(*->+„ and p_D<*)+ 7iBr 
are the CM energy and momentum, respectively, of the 
-D(*) + 7i sr combination. This peak is expected to be wide 
and asymmetric due to the photon energy resolution and 
higher-order corrections to ISR (i.e., more than one 7i sr 
in the event). From the Monte Carlo (MC) simulation the 
resolution of this peak is estimated to be ~ 300MeV/c 2 , 
which is not sufficient to separate the DD*, D*D* or 
£)(*) D*^tt final states. To disentangle the contribu- 
tions from these final states and to suppress combina- 
torial backgrounds, we use the slow pion from the unre- 
constructed D*~ . The difference between the mass re- 
coiling against D^*^ + 7i sr and -D < -*- l+ 7r s ^ ow 7i sr (recoil mass 



To measure the e + e~ 
smaller than the initial 

ergy {Eqm) at £?-factories, ISR can be used 
allows a measurement of cross sections in a broad energy 
range while the high luminosity of the _B-factories com- 
pensates for the suppression associated with the emission has a narrow distribution (a ~ 1.4MeV/c 2 ) around the 



ISR difference): 



AM rec =M rec (£>(*)+7 isr )-M rec (£>(*)+7rT 7isr ) , (2) 



3 



nominal m D *- — m-p , since the uncertainty in 7i sr mo- 
mentum partially cancels out. 

For the measurement of the exclusive cross section, 
one needs to determine the D*~ mass when one 

of the D*'s is not reconstructed. In the absence of 
higher-order QED processes, M(D^ + D*^) is the mass 
recoiling against the 7i sr . However, the photon energy 
resolution results in a typical M roc (7i sr ) resolution of 
~ 100 MeV, which is too wide for the study of rela- 
tively narrow 

D (*)+ D *- 

mass states. We significantly 
improve the M roc (7i sr ) resolution by applying a refit that 
constrains M rec (-D^*' )+ 7isr) to the nominal D*~ mass. 
In this way we use the well-measured properties of the 
fully reconstructed £>(*)+ to correct the energy of the 
7 isr . As a result, the M fl(<) ± c , T (= M rcc (7 isr )) res- 
olution is improved by a factor of ~ 10 and varies 
from ~ 6MeV/c 2 around threshold to ~ 12MeV/c 2 at 
M D (,) +D ,- = 5.0GeV/c 2 . The recoil mass difference af- 
ter the refit procedure (AM®.) has a resolution improved 
by a factor of ~ 2. Finally, the cross section is derived 
from the mass spectrum after the refit. 

All charged tracks should originate from the inter- 
action point (IP) with the requirements dr < 2 cm 
and dz < 4 cm, where dr and dz are the impact pa- 
rameters perpendicular to and along the beam direc- 
tion with respect to the IP. Charged kaons are required 
to have the ratio of particle identification likelihoods, 
V K = C K /{C K +Ar) 0, larger than 0.6. No identifi- 
cation requirements are applied for pion candidates. Kg 
candidates are reconstructed by combining 7r + 7r _ pairs 
with an invariant mass within 10 MeV/c 2 of the nomi- 
nal K s mass. The distance between the two pion tracks 
at the K s vertex must be less than 1 cm, the transverse 
flight distance from the interaction point is required to 
be greater than 0.1cm, and the angle between the K s 
momentum direction and the flight direction in the x — y 
plane should be smaller than 0.1 rad. Photons are recon- 
structed in the electromagnetic calorimeter as showers 
with energies greater than 50 MeV that are not associ- 
ated with charged tracks. ISR photon candidates are 
required to have energies greater than 2.5 GeV. Pairs 
of photons are combined to form ir° candidates. If the 
mass of a 77 pair lies within 15MeV/c 2 of the nominal 
7T mass, the pair is fit with a 7r° mass constraint and 
considered as a ir° candidate. 

D° candidates are reconstructed using five decay 
modes: K~ir + , K~K + , K~ir~n + ir + , Kgir + Tr~ and 
K~ir + ir°. A ±15MeV/c 2 mass window is used for all 
modes except for K~Tr~ir + ir + , where a ±10MeV/c 2 re- 
quirement is applied (~ 2.5 a in each case). D + candi- 
dates are reconstructed using the decay modes Kgir + , 
K~tt + tt + and K~K + tt + . A ±10MeV/c 2 mass window 
is used for all D + modes. To improve the momentum res- 
olution of D meson candidates, final tracks are fitted to a 
common vertex applying the nominal D° or D + mass as a 
constraint. D* candidates are selected via D* + — > D°tt + 



and D*° -> £>°7r° decay modes with a ±2 MeV/c 2 D* -D 
mass-difference window. A mass- and vertex-constrained 
fit is also applied to D* candidates. 

The distribution of M rcc (Z?* + 7i sr ) in the data, with- 
out any requirements on the slow pion, is shown in 
Fig. [T]a). The excess around the D*~ mass includes 
the e + e~ — > D* + D*~ y lsr signal as well as contributions 
from the e + e~ — * D + D*~y lsr channel. The shoulder 
at higher masses is due to e + e~ — > D* + D^*'TTji sr . The 
excess that is evident at ~ 2.5 GeV/ c 2 corresponds to 
e + e~ — > D* + D**~ji sr . The background from the other 
processes is substantially suppressed by the inclusion of 
the slow pion from the unreconstructed D*~ and the 
tight requirement on AM r *; i.e., within ±2MeV/c 2 of 
the nominal mj).- — m-^ , a clean peak corresponding 
to e+e~ -> D*+D*~ is evident in the M rcc (£>*+7 isr ) dis- 
tribution (Fig. Q]b)). We define the signal region by the 
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FIG. 1: The M rcc (D ( * )+ 7 isr ) distribution for the data: a), d) 
without requirement of slow pion; b), e) with AM* require- 
ment. Histograms show the normalized M D („)+ sidebands 
contributions, c), f) The distribution of AM r e' c in the data 
after the refit procedure. The selected signal windows are 
indicated by the vertical lines. 

requirement that M rec (D* + 7i sr ) be within ±0.2GeV/c 2 
of the nominal D *~ mass. This tight requirement sup- 
presses e + e~ — * D* + _D*~7r°7i sr events, which have a sim- 
ilar AM®, distribution. The spectrum of AM®, for the 
signal M rcc (Z?* + 7i sr ) window after the refit procedure in 
data is shown in Fig. [T]c). 

In case of multiple entries, the D* + ir candidate with 
the minimum value of Xtot — Xm(d°) ~^~Xm(d*+) +Xam b ' 



4 



is chosen, where each \ 2 is defined as the squared ratio 
of the deviation of the measured parameter from the ex- 
pected signal value and the corresponding resolution. 

The following sources of background are considered: 
(1) combinatorial background under the reconstructed 
D* + peak; (2) real D* + mesons coming from the signal or 
other processes combined with a combinatorial slow pion; 
(3) both the D* + and slow pion are combinatorial; (4) 
the reflection from the process e 



+ r,- 



D*+D 



''"missTisr 



with an extra ir° in the final state; (5) the contribution of 
e + e~ — * D* + D*~tt° when an energetic ir° is misidentificd 
as a single 7j sr . The contributions of backgrounds (1) 
and (2) are extracted using D*+ and AM*' sidebands, 
which have twice the area of the signal region. The latter 
sideband is shifted by 10MeV/c 2 from the signal region 
to avoid signal over-subtraction due to the higher-order 
ISR tail. Background (3) is present in both the M(D*+) 
and AM* sidebands and is, thus, subtracted twice. To 
account for this over-subtraction we use a 2-dimensional 
sideband region, where events are selected from both the 
D*+ and the AM*' sidebands. The total fraction of 
combinatorial backgrounds (1-3), found to be as large 
as ~ 10%, is subtracted from the signal-region D* + D*~ 
mass spectrum. 

Process (4) produces a broad peak in the AM*' distri- 
bution around the nominal m D ,- ~m-^ value and, thus, 
is not contained in the AM* sidebands. The dominant 
part of this background is suppressed by the tight require- 
ment on M rec (£>* + 7isr). The remaining part is estimated 
by applying a similar partial reconstruction method to 
the isospin-conjugate process e + e~ — > £>*°-D*~7r+ iss 7i sr : 
a D*° is fully reconstructed and the slow pion from a D*~ 
is used for the AM*' requirement. Since there is a charge 
imbalance in this final state, only events with a missing 
extra 7r+ iss can contribute to the AM r *' c peak. To extract 
the level of background (4) , this spectrum is rescaled ac- 
cording to the ratio of D*° and D* + reconstruction effi- 
ciencies and an isospin factor of 1/2. The contribution 
from background (4) is found to be consistent with zero. 
Uncertainties in this estimate are included in the system- 
atic error. The contribution of background (5) is deter- 
mined from reconstructed e + e~ — > D* + D*~ir° events us- 
ing a similar partial reconstruction technique but with an 
energetic 7r° replacing the 7; sr . The contribution of this 
background is found to be negligibly small; uncertainties 
in this estimate are also included in the systematic error. 

The analysis of the e + e _ — > D + D*~ cross sec- 
tion is identical to that described above for e + e~ — > 
£)*+£)*- w ith the fully reconstructed D* + meson re- 
placed by a fully reconstructed D + meson. The dis- 
tribution of M rcc (Z? + 7i sr ) with no requirements on the 
slow pion from the D*~ is shown in Fig. Q]d). The 
excess around the nominal D*~ mass corresponds to 
the e + e~ — * D + D*~y ml signal plus overlaps from the 
e + e~ — > D* + D*~ji SI channel. The shoulder at higher 
masses is due to e + e~ — > D*~ nji SI . The require- 



ment of a detected slow pion from the unreconstructed 
D*~ and a tight requirement on AM* provides the 
D + D*~ signal peak in the distribution 



clean e + e 



of M roc (D + 7i sl .) that is shown in Fig. Q]e). The AM*' C 
distribution for the signal M rec (D + 7i sr ) window is shown 
in Fig. [l]f). In the case of multiple entries, we apply a 
single candidate selection procedure similar to that used 
for D* + D*~ to all the distributions shown below. 

Similar sources of background (1-5) are considered for 
the e + e~ — > D + D*~ study. The contributions of back- 
grounds (1-3) are determined using D + and AM*' C side- 
bands with twice the area of the signal region. Here 
again, background (3) is present in both M(D + ) and 
AM*' sidebands and thus subtracted twice. To account 
for this over-subtraction we use a 2-dimensional sideband 
region that contains pure background (3) events. 

The level of contamination from background (4) is 
determined from isospin-conjugate events, e + e~ — > 
£>°£>*~7r+ iss 7 isr , in the data. The D + D*~ analysis is 
repeated with the fully reconstructed D +, s replaced by 
fully reconstructed D°'s. The D°D*~ mass distribution, 
after combinatorial background subtraction, is rescaled 
according to the ratio of D + and D° reconstruction ef- 
ficiencies and an isospin factor of 1/2. By chance, the 
contamination from the process e + e~ — > D* + D*~, fol- 
lowed by D* + — > D + tt q is also included in our estimate 
with correct scaling because B(D* + — > D + tt°)/B(D* + — > 
£>°7r + ) is also approximately 1/2. The contribution from 
background (5) determined from reconstructed e + e~ — > 
D + D*~tt° events in the data is found to be negligibly 
small and taken into account in the systematic error. The 
total background level is ~ 20% or less of the signal for 
all values of M(D+D*~). 

The e + e~ — > £)(*)+£)*- cross sections are extracted 
from the D*~ mass distributions by the rela- 

tion Em 



a(e H 



£)(*)+£)*-) = 



dN/dm 
rjtotdL/dm 



(3) 



where m = M(D^+D*-), dN/d m is the mass spec- 
tra obtained before corrections for resolution and higher- 
order radiation, while ?7tot is the total efficiency. The 
factor dL/dm is the differential ISR luminosity 



dL/dm = — ( (2 - 2x + x 2 ) In - 

TTX V 1 



1 + C 2 \ 2mC 

x O i 



c 



rp2 
°CM 



(4) 



where x = 1 — m 2 /Eq M , L is the total integrated luminos- 
ity and C — cos 0q , where 9q defines the polar angle range 
for 7 isr in the e + e~ CM frame: 9 a < 6> 7iar < 180 - 6 a . 
The reconstruction efficiencies, determined as a function 
of M(D^ + D*-) by MC simulation, are found to be in- 
dependent of M{D { ^+D*-) for both processes and are 
equal to r]{D*+D*-) = 4.3 x 10~ 3 and rj(D+D*-) = 
3.9 x 10~ 3 . The resulting exclusive e+e" D^+D*- 
cross sections are shown in Fig. [5] with statistical uncer- 
tainties only. Since the bin width is much larger than the 
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FIG. 2: The exclusive cross sections for a) e + e — > D* + D* 
and b) e + e~ -> D+D*-. The error bars include only the 
statistical uncertainties. 



resolution, no correction for resolution is applied. Since a 
reliable fit to the cross-sections obtained above requires a 
solution to a non-trivial and model-dependent problem of 
coupled channels and threshold effects, we do not report 
results here. 

The systematic errors for the er(e + e~ — * D^ + D*~) 
measurements are summarized in Table HI The system- 



TABLE I: Contributions to the systematic error on the cross 
sections, [%]. 



Source 


D*+D*~ 


D + D*- 


Background subtraction 


±5 


±4 


Angular distributions 


±5 




Reconstruction 


±7 


±6 


Cross section calculation 


±2 


±2 


B(£> w ) 


±3 


±5 


MC statistics 


±4 


±3 


Kaon identification 


±1 


±1 


Total 


±11 


±10 



atic errors associated with the background (1-3) subtrac- 
tion are estimated to be 3% from the uncertainty in the 
scaling factors for the sideband subtractions using fits to 
the M(DW + ) and AM^ distributions in the data with 
different signal and background parametrization. Un- 
certainties in backgrounds (4-5) are estimated conserva- 
tively to be smaller than 2% of the signal in the case of 
D* + D*~; these two sources are added linearly to give 5% 
in total. In the case of the D + D*~ , backgrounds (4-5) 
are subtracted using the data and only the uncertainty 



in the scaling factor for the subtracted distribution is 
taken into account. A second source of systematic error 
is due to the unknown hclicity angle composition of the 
D* + D*~ final state which can be a mixture of D^~D^T 
D^D*^ and D^D*^, where the subscripts L and T 
refer to longitudinally and transversely polarized D*'s. 
For the efficiency calculation, we assume equal fractions 
of these helicity states and consider the extreme cases 
(pure DljT and pure D*^~ D*^) for the efficiency un- 
certainty. There is no corresponding uncertainty in the 
case of the D + D*~ final state, where the D*~ polariza- 
tion is fixed by angular momentum and parity conserva- 
tion. A third source of systematic error comes from the 
uncertainties in track and photon reconstruction efficien- 
cies, which are 1% per track, 2% per slow pion and 1.5% 
per photon, respectively. The systematic error ascribed 
to the cross section calculation is estimated from a study 
of the C-dependence of the final result and includes a 
1.5% error on the total luminosity. Other contributions 
come from the uncertainty in the identification effici ency 
and the absolute D° and £)(*)+ branching fractions [14| . 

In summary, we report the first measurements of ex- 
clusive e + e~ — > D* + D*~ and e+e~ — > D + D*~ cross 
sections at y/s around the D* + D*~ and D + D*~ thresh- 
olds with initial state radiation. The shape of the 
e + e~ — > D* + D*~ cross section is complicated with sev- 
eral local maxima and minima. The minimum near 
4.25GeV/c 2 — in the Y(4260) region — could be due to 
D*D* (DD**) threshold effects described in [IE EE Ell 
or due to destructive interference of this state with other 
ip(nS) states. Aside from a prominent excess near the 
-0(4040), the e + e~ — * D + D*~ cross section is relatively 
featureless. The measured cross sections are compati- 
ble [HI within errors with the D^D* exclusive cross 
section in the energy region up to 4.260 GeV measured 
by CLEO-c 0. 
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